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ABSTRACT 

Recent observations reveal galaxies in the early Universe (2 < z < 6.4) with large reservoirs of 
molecular gas and extreme star formation rates. For a very large range of sources, a tight relationship 
exists between star formation rate and the luminosity of the HCN J=l-0 spectral line, but sources at 
rcdshifts of z^2 and beyond do not follow this trend. The deficit in HCN is conventionally explained 
by an excess of infrared (IR) radiation due to active galactic nuclei (AGN). We show in this letter 
not only that the presence of AGN cannot account for the excess of IR over molecular luminosity, but 
also that the observed abundance of HCN is in fact consistent with a population of stars forming from 
near-primordial gas. 

Subject headings: stars: formation — galaxies: abundances — galaxies: starburst — astrochemistry 



1. INTRODUCTION 

Recent work (Gao & Solomon 2004, Wu et al. 2005) 
indicates that there exists a strong correlation between 
the luminosity of molecular lines - particularly the HCN 
J=l-0 transition- and the star formation rate derived 
from infrared luminosity. Impressively, these relations 
apply to individual star forming regions in the Milky Way 
and also to external systems across ten orders of magni- 
tude in luminosity. However, the three highest redshift 
systems in the sample used by Gao & Solomon (at z ~ 2) 
do not fit this trend, but have luminosities in HCN lower 
than would be expected given their star formation rate. 
A detection of the HCN J=5-4 line in APM 08279+5255 
at a redshift z=3.91 has been reported with Lpjn/ Lhcn 
similar to that seen in other high-z sources (Wagg et al. 
2005). Limits on HCN luminosities obtained for systems 
at much higher redshifts (z ~ 6) (Carilli et al. 2002) 
indicate that these systems also do not fit the low-z re- 
lation. 

The conventional explanation for these discrepancies 
is to attribute the observed infrared excess to AGN ac- 
tivity. The z ~ 2 systems are, indeed, quasar hosts, 
presenting us with the problem of distinguishing IR ra- 
diation associated with star formation and that due to 
the AGN. Detailed modeling of the spectral energy dis- 
tribution (SED) of one such source, the Cloverleaf quasar 
(Weifi et al. 2005), has shown that there exists a warm 
dust component which is believed to be associated with 
star-forming activity. In this letter we test the hypothesis 
that hidden AGN activity is responsible for the deviation 
of these high-z sources from the observed relationship be- 
tween infrared and HCN luminosity. By using chemical 
models of high-mass star formation, we in fact find that 
AGN activity should be accompanied by an increase in 
the relative abundance of HCN. Furthermore, we find 
that models of star formation from near-primordial gas, 
enriched only by the ejecta from population III stars, 
can account for the observed HCN abundance; hence we 
show that there is no 'lack' of HCN at high redshift. 
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2. THE EFFECT OF AGN 

Our previous work (Lintott et al. 2005) used a time- 
dependent chemical model to simulate the formation of 
massive stars in order to constrain the expected emission 
from molecular species at high redshifts. Here, we make 
use of the same model to explore the effect of AGN on 
high redshift star forming systems. We use a model with 
initial atomic abundance ratios from calculations of the 
yield from zero-metallicity supernovae (Chieffi & Limongi 
2004) and overall metallicity l/100th solar. We then 
model the initial collapse of the protostellar clump, in- 
cluding chemical reactions both in the gas and on the sur- 
faces of dust grains. Once a critical density is reached, we 
assume the high-mass protostar ignites, causing species 
trapped in the icy grain mantles to sublimate. The newly 
released molecules radiate strongly in the warm gas; the 
protostar has become a 'hot core'. Although short-lived, 
hot cores are among the most luminous of all sources of 
sub-mm radiation in the Milky Way. Following previ- 
ous work on low-mass star formation (Lepp & Dalgarno 
1996) we assume that the effect of AGN activity within 
the host galaxy will be an increase in ionization flux due 
to x-rays. 

Ionization is of crucial importance in driving the chem- 
istry of star formation regions; the most significant gas 
phase reactions within dark clouds, particularly for oxy- 
gen and nitrogen bearing species, are ion-neutral reac- 
tions. We simulate the increase in x-ray ionization by 
enhancing the standard (1.0 x 10 17 s _1 ) cosmic ray ion- 
ization rate by a factor, (, of 2 to 1000. The results of 
our modeling are shown in Figure 1. 

The primary formation route for HCN in the bench- 
mark model is via the reaction of H2 and CN. The lat- 
ter is formed from atomic nitrogen and C2. UV pho- 
tons formed as a by-product of the increased ionization 
produce enhanced abundances of atomic nitrogen, pro- 
moting the formation of CN and hence HCN. In addi- 
tion, the increased ionization promotes the formation of 
HCN from dissociative recombination of HCNH+ ; this is 
the dominant formation route for HCN during the early 
stages of models with increased ionization. Note that 
reactions on the surface do not affect the abundance of 
HCN in our model as there is no significant formation 
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Fig. 1. — Change in abundance for HCN relative to a 'bench- 
mark' model with only cosmic ray ionization and Milky Way metal- 
licity and atomic abundances. The solid line represents a model 
with only cosmic ray ionization and abundances from Chicffi and 
Limongi (2004), standardized to l/100th solar. Time since the ig- 
nition of the protostar and consequent sublimation of species from 
grain mantles shown on the x-axis. 



or destruction occuring on the grains. However, if freeze 
out is less efficient, more HCN would form in the gas 
phase during collapse. 

A £ of 100 results in an increase in relative abundance 
in HCN at all but very early times. Runs with £ = 1000 
produce an initially higher fractional abundance than 
those with ( — 100, although the final abundance is un- 
affected. This result echoes that of Lepp and Dalgarno in 
low-mass star forming regions, although they find that, 
in these sources, the abundance of HCN (in common with 
all species in their model) declines at the high ionization 
rates explored by our model. In the model described here 
the increased density during collapse and the resulting 
freeze-out of molecules prevents dissociation. 

Our models show clearly that an increase in HCN 
abundance is likely to accompany AGN activity. The 
magnitude of the increase (for relatively modest increases 
in ionization rate) can be as large as two orders of magni- 
tude. The far-IR excess seen in the most distant sources 
included by Gao & Solomon is approximately an order 
of magnitude (a typical source with log(LncN)^ 9.5 has 
log {Lib) ~ 13.5 while their relation would have pre- 
dicted log(Ljfl) ~ 12.5). Note that we do not model 
the physics of the AGN but only its probable chemical 
effects, and further research would be necessary to in- 
vestigate the relationship between the increases in ion- 
ization and infrared luminosity. However, it seems likely 
that the consequent increase in HCN abundance from a 
AGN powerful enough to cause an order of magnitude 
increase in IR luminosity is large enough that the devi- 
ation from the observed trend cannot be explained by 
the presence of AGN. Similar enhancements of HCN in 
nearby AGN are observed (Kohno et al. 2001), and the 
use of HCN/HCO++ and HCN/CO ratios as diagnostics 
for AGN activity has been suggested. In addition, our 
preliminary results suggest that other species may also 
be worth investigating, in particular CS and SO2. 

3. IR EXCESS OR NITROGEN UNDER-ABUNDANCE? 

If not AGN activity, what could be responsible for the 
deviation of molecular luminosities in high redshift sys- 
tems from a relation which holds over a large range of 



luminosities? In a previous paper (Lintott et al. 2005) we 
suggested that the molecular gas in systems with rapid 
rates of star formation could be in the form of hot cores. 
Furthermore, gas enriched only by large zero-metallicity 
stars should be under-abundant in nitrogen compared to 
the solar neighbourhood. Hot cores formed from such gas 
should therefore have lower than expected abundances of 
nitrogen-bearing compounds. 

In this section, we investigate whether the molecular 
abundances derived in Section 2 are consistent with ex- 
isting detections of HCN in high-redshift galaxies, and 
thus the hypothesis that deviations from the standard 
LyrJ -^hcn ratio are due not to an excess of infrared lu- 
minosity, but an under-abundance of HCN due to star 
formation from near-primordial gas. 

We consider two sources observed using the Very Large 
Array (VLA), the Clovcrleaf quasar (H1413, z=2.56) 
(Solomon et al. 2003) and J1409 (z=2.58, Carilli et al. 
2002). Unfortunately, both of these sources have signif- 
icant optical depths. However, we can still derive the 
column density as long as we assume a beam filling fac- 
tor and excitation temperature (Nguyen-Q-Ricu et al. 
1992). As the SED of the Cloverleaf has been studied in 
detail and a warm component with Td = hQK associated 
with star formation detected, we will assume T ex = 50K 
for both sources. We estimate the beam filling factor 
by assuming that the only source of HCN emission arises 
from hot cores with radii of 0.3pc. We have assumed such 
objects are larger than their Milky Way analogues in or- 
der to maintain an extinction large enough to harbour 
substantial reservoirs of molecular gas. 

Our model predicts the contribution of each individ- 
ual core to the column density, and, thus, given the 
number of cores in the source the overall column den- 
sity. The observed luminosities for the Cloverleaf and 
J1409 in K kms" 1 pc 2 are L HCN = 3.2 ± 0.5 x 10 9 and 
Lhcn = 6.7 ± 2.2 x 10 9 respectively. 

We find that our model is consistent with the obser- 
vations in both cases if approximately 10 7 hot cores are 
present. The column densities derived for the two sources 
are then 1.8 x 10 17 cm~ 2 (Cloverleaf) and 3.8 x 10 16 cm~ 2 
(J1409). In Lintott et al. 2005 we estimated that the 
Milky Way contains 10 5 hot cores and assert that the 
number of hot cores scales with star formation rate. As 
both of these sources have star formation rates (as de- 
termined by the infrared luminosities) greater than two 
orders of magnitude above the Milky Way rate, 10 7 hot 
cores seems realistic. It should be noted that we have 
assumed that the observed molecular gas is entirely in 
the form of hot cores, and that these values should be 
considered as upper limits. It is clear from Figure 1 that 
the relative abundance of HCN in these models is at least 
a factor of 10 4 below that expected in hot cores formed 
from gas with solar abundance ratios. Even if the num- 
ber of hot cores in these rapidly star-forming systems is 
the same as in the Milky Way, hot cores with solar abun- 
dance ratios would result in a luminosity of HCN greatly 
in excess of the observed value. 

Thus there is no discrepancy between the star forma- 
tion rate (derived from the IR luminosity) and the ob- 
served luminosity of HCN, if stars are forming from near- 
primordial gas at redshifts of z=2-4. 

4. CONCLUSIONS 



Rapid star formation in the presence of active galactic nuclei 
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We have argued that the observed differences between 
£ir/£hcn at low and high rcdshifts cannot be due to the 
presence of AGN enhancement of the infrared luminosity. 
Such activity, we show, would increase the abundance of 
HCN in the observed systems. Detailed chemical models 
of such systems show promise in distinguishing between 
AGN enhanced and purely star forming activity through 
the ratio of HCN with species which are not enhanced 
by AGN. Examples of such species may include CS and 
SO2 (both of which are enhanced in galactic hot cores 
(Hatchell et al. 1998) and are expected to be favoured 
in nitrogen-deficient gas (Lintott et al. 2005). 

We explain the observations instead by hypothesiz- 
ing that stars are forming from near-primordial mate- 
rial. Such material is believed to be deficient in nitrogen 
compared to the solar neighbourhood. 

Jirnincz and Haiman (2006) discuss four observational 
constraints on the nature of star formation at redshifts 
between 3 and 4. Specifically, they explain significant 
UV emission in Ly-break galaxies, Ly-alpha emission iso- 
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lated from ionizing continuum radiation, a population of 
sources with strong Ly-alpha emission and strong Hell 
lines by proposing a significant fraction of star formation 
from primordial gas. These results imply that mixing of 
metals is inefficient within galaxies, and so solar abun- 
dance ratios may not have been reached at a redshift of 
z=2. 

Such models present a challenge to existing models of 
chemical enrichment, which assume that the first popu- 
lation II stars form at much higher redshifts. Simulations 
(e.g. Abel et al. 2002) which suggest that population III 
stars will be massive assume that such stars will form at 
rcdshifts much greater than 4; it may be that the predic- 
tions of yields from such stars may need to be revised. 
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Reaction 



H + HNC — ► HCN + H 
N + CH 2 — ► HCN + H 

H 2 + CN — ► HCN + H 
N + HCO — ► HCN + O 
HCNH++e — ► HCN + H 

H + H 2 CN — > HCN + H 2 



H+ 



HCN — > HCNH+ + H 2 
HCN — ► HCN+ + H 
H3O+ + HCN — > HCNH+ + H 2 



HCN + Photon 



CN + H 



Rate coefficient /cm 3 s ~ 



..2 X 10" 10 exp — 



300K , 
, 2.87 



4.04 x 10-13 (_£_) ' cxp - 
1.70 x 10" 10 
9.0xl0- 8 (3^) _0 ' 5 



1.0 x 10" 



' T V 
300K I 



8.1 x 10" 9 
1.05 x 10- 8 (3^ 
4.0 x 10 -9 



-0.13 



1.3 X lQ- 9 exp(-2.1A„) 



TABLE 1 

Major formation and destruction reactions of HCN; all those which account for 10% of the formation or destruction 

RATE AT ANY POINT IN BOTH BENCHMARK AND ENHANCED IONIZATION MODELS HAVE BEEN INCLUDED. THOSE IN BOLD ARE IMPORTANT ONLY 
IN MODELS WITH ENHANCED IONIZATION (ALTHOUGH NOTE THE FINAL REACTION IS IMPORTANT ONLY IN THE VERY EARLY STAGES). RATES 
ARE TAKEN FROM THE UMIST99 RATE FILE. T IS THE TEMPERATURE AND A v IS THE EXTINCTION. 



